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1.  Introduction 

This  report  describes  a  mathematical  model  to  calculate  aural 
detection  of  low  flying  aircraft.  Instead  of  trying  to  modify  an 
existing  model  computer  program  [1],  it  was  decided  to  construct  a  new 
program  based  on  an  8-step  method  [2,3).  Guidelines  were  to  develop  a 
model  with  one  observer  and  one  aircraft. 

The  main  purpose  of  the  new  model  program  (Figure  1)  is  to  compute 
the  probability  of  aural  detection  with  time.  Aural  detection  is  based 
on  an  incremental  change  in  loudness  level  above  the  ambient  level.  The 
loudness  level  is  a  method  of  analysis  for  complex  sound  by  relative 
judgments  to  a  1000-nertz  tone-.  The  units  of  loudness  and  loudness 
level  are  sones  and  phons,  respectively,  and  will  be  described  further 
in  Section  3. 

One  item  of  concern  in  the  model  is  in  the  human  detection  portion. 
The  data  used  to  simulate  human  detection  were  based  upon  two  subjects. 
If  available,  better  data  should  be  incorporated  into  the  model  to  give 
a  truer  reoresentation  of  a  human  observer. 


2.  Prediction  Procedure 

In  predicting  the  noise  level  which  will  be  received  by  the 
observer  on  the  ground  from  an  airborne  source,  it  is  necessary  to  know 
the  following: 

a)  The  sound  pressure  level  (SPL)  at  a  known  reference 
distance  from  the  source 

b)  The  distance  from  the  source  to  the  point  of  observation 
on  the  ground 

c)  The  sound -propagation  characteristics  cf  the  atmosphere 
along  the  path  between  the  source  and  the  ground  point  [2] 

When  measuring  the  sound  signature  of  the  aircraft  at  a  known 
distance,  the  distance  should  be  close  to  the  aircraft  to  limit 
atmospheric  effects  on  the  propagation  of  the  sound.  Aircraft  sounds 
are  directional  and  must  be  treated  as  such  when  taking  measurements  at 
a  known  reference  distance  [2?.  The  SPL  at  the  observer  is  calculated 
by: 


SPL  =  SPL  -  20  log 


10  (“D  ' 


dB  re  0.0002  dyne -'em” 


1 


where 


SPL^  =  reference  sound  pressure  level  (dB) 

d  =  distance  from  source  at  which  reference  level  is 
measured 

d  s  distance  from  source  to  observer 

A  =  total  attenuation  (d3)  over  distance  a  -  d  . 

o 

The  term,  20  log-^d/d^) ,  in  Equation  (1)  is  the  reduction  in  SPL 

due  to  spherical  divergence  of  the  sound  wave  in  spreading  outward  from 
the  source,  commonly  referred  to  as  "inverse  square  law  attenuation." 

If.  is  equivalent  to  a  6-decibel  reduction  in  SPL  for  each  successive 
doubling  of  distance  from  the  source  £2].. 

Total  attenuation.  A,  is  determined  by  three  types  of  attenuation 
constants:  (a)  atmospheric  absorption  due  to  humidity,  (b)  atmoshperic 
abosrption  in  excess  of  humidity,  in  particular  scattering  of  sound 
waves  due  to  inhomogeneity  and  turbulence  in  the  atmosphere,  and 
(c)  attenuation  due  to  terrain  absorption. 

A  -  *h  +  Ao  •  (2) 

where 

Ag  =  atmospheric  absorption  due  to  humidity 

A  -  residual  attenuation  (atmospheric  absorption  in 
excess  of  humidity  attenuation) 

=  terrain  absorption. 

Humidity  attenuation,  A_,  can  be  calculated  by  the  formula 

n 

expressed  as  an  attenuation  coefficient,  a„,  in  decibels/1000  feet  f 2] : 

ri 


^  =  O.X  f  (T  -r  45)/((fAn“)  +  (hVf)] 

d,dQ  =  slant  distance  (1000  feet) 

f  =  geriatric  mean  frequency  of  an  octave  band  (kilocycles)  , 

T  =  temperature  (=F). 

Charts  of  plotted  against  absolute  humidity  and  temperature  for  various 
octave  frequency  bands  are  given  in  Figure  2  [2} . 

Aq,  residual  attenuation,  is  due  principally  to  the  scattering 

affects  of  air  turbulence.  Calculation  of  Aq  depends  on  whether  the 

observer  is  in  a  sound  shadow.  The  shadow  region  is  due  to  upward 
refracted  sound  rays  and  corresponds  to  a  region  in  which  sound  is 
highly  attenuated.  This  refraction  is  due  to  a  decrease  in  sound 
velocity  with  increasing  altitude.  To  calculate  whether  shadow  forma¬ 
tion  has  occurred,  the  following  procedure  is  followed.  The  sound 
velocity  in  a  given  horizontal  direction  in  the  atmosphere  depends  on 
the  temperature  and  on  the  vector  addition  of  the  cosponent  of  wind 
velocity  in  that  direction.  The  sound  velocity,  c,  is  given  by  the 
expression: 


c  =  c„  -T-  c 
i  K 


wnere 


the  sound  velocity  determined  by  temperature  at 
the  given  ooint 


%  =  ^30  cos  <6> 


=  the  wind  speed  at  30  feet  height 
C  =  the  wind  sound  angle  (Appendix  A) 


W„  = 
30 


1  :  0.2?  iOgiQ 


?_  =  the  known  wind  magnitude  at  height,  z. 


Sound  propagation  in  directions  for  which  ct,  is  positive  is  referred  to 

as  "downwind"  propagation  and  "upwind"  propagation  for  negative  c,^. 

Khan  the  variation  of  cT  and  with  height  is  such  that  the  net  sound 

velocity,  c,  decreases  with  height,  sound  rays  will  be  refracted  upward 
as  shown  in  Figure  3>.  To  determine  if  the  observer  is  in  a  sound 
shadow,  use  Figure  4  to  find  R^,  the  horizontal  distance  from  source  to 

shadow  boundary.  This  calls  for  the  calculation  of  the  logarithmic 
gradient,  B,  defined  by: 


3  =  Bt  -r  a. 


a,  =  4Kx  l<f  4 
2, 


10810  (t) 


where 


T,Tq  =  temperature  at  target  and  observer 


~-,~o  -  height  c-f  target  and  observe 


=  1.5  cos  (e>  10 


If  the  horizontal  distance  from  source  to  receiver,  r,  is  greater  than 

3  R  ,  the  observation  point  will  be  in  a  sound  shadow.  The  calculation 
0 

of  A  when  the  observer  is  in  a  sound  shadow  fellows.  For  wind  sound 
o 

angles  between  110  and  150  degrees,  use  a  value  of  30  decibels  for 
residual  attenuation,  A  -  Otherwise,  determine  A  at  the  distance, 

R^,  by  Figure  5  (the  curves  have  been  extrapolated  by  the  user  for  more 

versatile  range  requirements)  and  add  an  amount  delta  Aq,  given  by; 


%  -  [(r  -  so)(30  -  V}]  2 


(ID 


■ 


—  V'' 


Average  values  of  residual  attenuation,  Aq,  due  to  turbulence  outside  a 

sound  shadow  are  plotted  in  Figure  5.  Weather  variables  affect  the 
estimation  of  residual  attenuation  for  source  elevation  angles  less  than 
2  degrees  for  all  frequency  bands  and  for  elevation  angles  between  2 
and  10  degrees  for  the  75- to  150-hertz  band.  Therefore,  under  these 
conditions  average  residual  attenuation  muse  be  corrected  by  the 
formula : 


wnera 


A.  *  A  -  0.3  b_  OVT  -5-  4.5)  -f  0.12  h.  (H, 

O  O  A  w  - 


Ao  =  predicted  value  of  attenuation 

Aq  =  average  value  given  by  curves  of  Figure  5. 

4*T  =  temperature  gradient  (ep/1000  feet) 

3T  =  2  ;t5  -  T  ; 


where 

T_  »  known  temperature  at  5  feat  height 
-  temperature  at  500  feet. 

The  coefficients,  and  h.,.  are  given  by  Table  1.  The  temperature 

a.  2: ' 

gradient,  Aa,  is  taken  as  positive  when  temperature  at  5  feet  is  higher 
than  that  at  500  feet  {2; . 

Aj.  allows  the  introduction  of  four  types  of  terrain  conditions  for 

simulation:  (1)  open  area.  (2)  18- inch  grass,  {3)  lightly  treed,  and 
(4)  forest.  Terrain  absorption  is  dependent  upon  elevation  angle, 
range,  agd  frequency.  The  farther  devnrange  and  the  smaller  the 
elevation  angle,  the  higher  the  absorption.  Figure  6  and  Tabic  2  show 
percent  of  terrain  absorption  coefficient  versus  elevation  angle  where 
100-psrcent  terrain  absorption  attenuation  corresponds  to  0-aegree 
elevation  angle.  To  utilise  Figure  6  end  Table  2,  read  a  percentage 
factor  from  Figure  6  then  multiply  this  factor  by  the  appropriate 
terrain  coefficient  of  Table  2.  The  corrected  coefficient  will  be  in 
terms  of  tHei&els/lCOO  feet  ;  1} - 


3.  Hearing  and  Detection 


Noises  differ  in  intensity  at  different  frequencies.  In  aural 
detection  one  must  consider:  (a)  ambient  noise  level,  (b)  signal  embedded 
in  the  ambient  level,  arid  (c)  the  sensitivity  of  the  ear  to  (a)- and  (b). 
One  must  compute  the  probability  of  detecting  an  incremental  intensity 
level  as  the  signal  comes  into  the  ambient  level.  To  approach  the 
problem  in  the  above  fashion,  loudness  levels  have  to  be  calculated; 

(c)  loudness  level  of  ambient  alone  and  (b)  loudness  level  when  a 
signal  has  been  added  to  the  ambient  level.  Loudness  levels  pertain  to 
the  magnitude  of  the  auditory  sensation  that  a  person  experiences  and 
should  not  be  confused  vith  other  distinguishing  characteristics  ot 
sound  (pitch,  penetrating  or  dull,  etc.).  Loudness  and  loudness  levels 
are  distinguished  and  measured  in  terms  of  sones  and  phons,  respectively. 
Acousticians  have  chosen  a  1000-hertz  tone  for  the  standard  reference  of 
the  phon  and  sone.  The  phon  level  of  an  unknown  sound  is  the  SPL  of 
a  1090-hertz,  tone  that  has  been  judged  by  listeners  to  equal  the 
loudness  of  the  unknown  sound.  The  sone  is  the  loudness  adjudged  to  a 
40-decibel,  1000-hertz  tone  [4] .  One  sone  equals  40  phons  (Fig:* 
and  is  related  by  the  equation: 


log. n  N  =  0.03  L  -1.2  (13) 

iu  n 

where  N  is  the  lojdnesu  (sones)  and  L^  is  the  loudness  level  (phons). 

With  t  le  method  of  calculating  loudness  in  sones  available,  the  stage 
is  set  to  calculate  probability  cf  detection: 

a)  To  calculate  loudness  of  background  alone  (sones): 

1)  Calculate  SPL  in  each  octave  band 

2)  Convert  to  sones  in  each  octave  band 

3)  Multiply  all  except  the  maximum  sone  value 
found  by  0.3 

4)  Convert  to  phons  by  Equation  (13)  (this  will  put 
the  loudness  level  in  terras  of  decibels  of.  a  1000- 
hertz  tone) . 

b)  To  calculate  loudness  of  signal  plus  ambient: 

1)  Calculate  SPL  of  signal  plus  ambient  in  each 

octave  band  (this  involves  changing  both  signal 
and  ambient  to  microbnrs,  sunaaing,  and  converting 
the  sum  back  to  decibels) 

Convert  to  sones  in  each  band  using  Figure  8. 


2) 


3)  Multiply  i.ll  except  the  maximum  sone  value 
found  by  0.3 

4)  Convert  to  phons  by  Equation  (13). 

These  calculations  will  give  decibel  intensity  values  of  sounds  in 
tef":s  of  a  1000-herts  tone.  Thus,  the  delta  decibel  value  above  the 
ambient  loudness  level  can  be  calculated  (Table  3) .  Using  the  Quantal 
Theory  of  Discrimination  for  detection  of  an  incremental  decibel  when 
the  increments  are  added  to  a  continuous  stimulus  [5J,  it  has  been  found 
that  the  listener  finds  it  difficult  to  distinguish  one-quantun  changes 
in  the  stimulus  from  the  changes  which  are  constantly  occurring.  The 
difference  reported  50  percent  of  the  time  is  equivalent  to  1.5  times 
the  quantal  increment  (decibels).  A  difference  reported  100  percent  of 
the  time  corresponds  to  2.0  times  the  quantum  increment  [5].  Table  4 
shows  the  mean  quantum  increment  for  two  persons  in  dec'bcls  as  a  function 
cf  sensation  level.  The  sensation  level  used  is  the  ambient  loudness 
level.  A  normal  distribution  is  used  to  find  a  probability  of  detecting 
a  delta  decibel  increment.  It  is  divided  into  12  regions  using  the 
mean  and  standard  deviations  of  Table  4.  These  data  were  plotted  and 
smoothed  out  to  give  a  more  uniform  curve.  The  standard  deviations  ere 
based  on  the  smoothed  data.  The  regions  correspond  to  the  proportion 
of  the  listening  population  assumed  to  have  that  quanta)  increment. 

A  weighting  factor  (Table  5)  describes  a  probability  of  an  individual 
being  in  tha  ~egion.  Using  the  mean  and  standard  deviations  of  Table 
4,  calculate  che  quantum  increment  of  the  midpoint  of  each  of  the  12 
regions.  To  calculate  the  probability  P,  let 


R’  = 


1 

R 

9 


if  R  is  less  than  or  equal  to  1 

if  R  is  between  1  and  2 

if  R  is  greater  than  or  equal  to  2 


then 


P  =  R'  -  1 


Starting  at  the  left  of  the  distribution  multiply  each  probability  by 
the  corresponding  weighting  factor  given  in  Table  5,  and  sura  up; 

12 

Probability  ^  =  Z  fv.  x  weighting  factor,  i) 

" ' '  ‘  i=l  '  f 


--**s 


4.  Summary 


In  predicting  the  SPL  that  will  be  received  by  the  observer  on 
the  ground  due  to  an  airborne  source,  several  factors  must  be  considered. 
Step  1  is  the  target  signature  measured  at  a  known  reference  distance 
at  aspect  angles  around  the  aircraft.  Conditions  which  must  be  known 
include:  (a)  temperature  at  observer  along  with  a  temperature  gradient, 
(b)  wind  velocity  vector  at  a  known  height,  (c)  the  type  of  terrain 
(there  are  four  types  available  for  simulation  in  the  model),  (d)  relative 
humidity,  and  (e)  the  target  flight  path.  To  estimate  the  SPL  received 
by  the  observer,  attenuation  along  the  acoustic  path  must  be  accounted 
for.  There  are  three  types,  of  attenuations  in  the  model:  (a)  attenua¬ 
tion  due  to  atmospheric  absorption  due  to  humidity,  (b)  attenuation  due 
to  terrain  absorption,  and  (c)  attenuation  due  _o  atmospheric  absorption 
in  excess  of  humidity,  usually  referred  to  as  residual  attenuation.  In 
determining  the  value  of  residual  attenuation  one  must  consider  if  the 
observer  is  in  a  sound  shadow.  Shadow  formation  is  due  to  upward 
refraction  of  sound  rays,  causing  a  highly  attenuated  region  where 
detection  is  poor  at  bast.  The  refraction  is  due  to  wind  and  temperature 
effects , 

At  this  point  the  estimated  SPL' s  of  the  ambient  and  signal 
received  by  the  observer  are  known  in  the  octave  bands.  The  number  of 
decibels  of  ambient  plus  signal  is  combined  in  each  octave  band  by 
conversion  of  microbars,  suraning  the  two,  and  converting  the  sum  back  to 
decibels.  The  loudness  in  sones  of  ambient  and  ambient  plus  signal  can 
be  obtained  at  this  stage  by  a  simple  mathematical  procedure  and  the  use 
of  a  table.  The  transformation  from  sones  to  phons  is  obtained  by  an 
equation.  Since  phons  are  equal  to  decibels  of  a  1000-hertz  tone,  the 
delta  decibel  can  be  calculated.  Using  t*  delta  decibel  and  a  derived 
normal  distribute  n,  a  probability  of  aurai  etection  is  obtained. 
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TABLE  3.  EXAMPLE  OF  THE  PREDICTION  OF  LOUDNESS  LEVEL* 


TABIE  5.  LIGHTING  FACTORS 


Example: 

Sensation  level  =  82  dB 

1. 

Mean  =0.285 

Standard 

2. 

Deviation  =  0.035 

3 

->  # 

Calculated 

delta  dB  =  0.27 

4. 

/» 

3. 

Region 

No. 


Midpoint 


0.1875 

0.20625 

0.22375 

0.24125 

0.25875 

0.27625 

0.29375 

0.31125 

0.32875 

0.34625 

0.35375 

0.38125 


standard  deviations  (12  regions) 

2.  Assign  weighting  factors 

3.  Calculate  midpoints  of  the 
intervals 

4.  Set  up  ratios 

5.  Check  to  see  that  the  ratios  do 
not  exceed  2.0  or  are  below  1.0 

6.  Subtract  1.0  from  each  ratio 
to  get  P. 

7.  Multiply  by  weighting  factors 

8.  Sum  up 


Quantals 

Ratio 

(R) 

?i 

Weighting 
Factor  (W /) 

Weighting 
Factor  x 

1.44 

0.44 

0.003 

0.0022 

1.31 

0.31 

0.017 

0.00527 

1.21 

0.21 

0.044 

0.00924 

1.12 

0.12 

0.092 

0.01104 

1.05 

0.05 

0.150 

0. 0007 5 

mSm 

0.00 

0.192 

0.0 

0.00 

0.192 

0.0 

JIB 

0.00 

0.150 

0.0 

0.821 

0.00 

0.092 

0.0 

0.779 

0.00 

0.044 

0.0 

0.742 

0.00 

0.017 

0.0 

0.708 

0.00 

.0.005 

0.0 

i 

Z  <W.  X  Pi) 

I 

=  0.02886* 

Probability  of  detection 


TABIE  5.  CONTINUED 


Example : 

Sensation  level  =  35  dB 

Mean  =  0.295 

Standard 

Deviation  =  0.025 

Calculated 

delta  dB  =  0.46 

i _ 

1.  Set  up  distribution  in  one-half 
standard  deviations  (12  regions) 

2.  Assign  weighting  factors 

3.  Calculate  midpoints  of  the 
intervals 

4.  Set  up  ratios 

5.  Check  to  see  that  the  ratios  do 
not  exceed  2.0  or  are  below  1.0 

6.  Subtract  1.0  from  each  ratio 
to  get  P. 

7.  Multiply  by  weighting  factors 

8 .  Sum  up 

Region 

No. 

Region 

Midpoint 

Quantals 

Ratio 

GO 

Pi 

Weighting 
Factor  (W.) 

Weighting 
lector  x  P^ 

1 

0.22625 

2.033149 

1.0 

0.005 

0.005 

2 

0.23875 

1,92670 

0.9267 

0.017 

0.01575 

3 

0.25125 

i:8308 

0.8308 

0.044 

0.0365 

4 

0.26375 

1.744 

0.7440 

0.092 

0.6844 

5 

0.27625 

1.665 

0.6650 

0.150 

0.09975 

6 

0.28875 

1.593 

0.5930 

0.192 

0.113856 

7 

0.30125 

1.5269 

0.5269 

0.192 

0.1011648 

8 

0.31375 

1.46613 

0.4661 

0.150 

0.069915 

9 

0.32625 

1.4099 

0.4099 

0.092 

0.0377108 

10 

0.33875 

1.3579 

0.3579 

0.044 

0.0157476 

11 

0.35125 

1.309 

0.3090 

0.017 

0.005253 

12 

0.36375 

1.2646 

0.2646 

0.005 

£  <»*  x  pjL) 

_ 

0.001325 

=  0.570410* 

_ 

Probability  of  detection. 


TABLE  5 .  CONCLUDED 


Example: 

Sensation  level  =  20  dB 

Mean  =  0.025 

Standard 

deviations  =  0.025 

Calculated 

delta  d3  =■  0.50 

1.  Set  up  distribution  is  one-half 
standard  deviations  (12  regions) 

2.  Assign  weighting  factors 

3.  Calculate  midpoints  of  the 
intervals 

4.  Set  up  ratios 

5.  Check  to  see  that  the  ratios  do 
not  exceed  2.0  or  are  below  1.0 

6.  Subtract  1.0  from  each  ratio 
to  get  P. 

7.  Multiply  by  weighting  factors 

8 .  Sum  up 

Region 

No.- 

Region 

Midpoint 

Quantals 

Ratio 

(R) 

P. 

l 

Weighting 
Factor  (W.) 

~  i 

Weighting 
Factor  x 

1 

0.29625 

1.6877 

0.6877 

0.005 

0.0034385 

2 

0.30875 

1.6194 

0.6194 

0.017 

0.0105298 

3 

0.32125 

1.5564 

0.5564' 

0.044 

0.0244816 

A 

0.33375 

1.498 

0.4980 

0.092 

0.459816' 

5 

0.34625 

1.444 

0.4440 

0.150 

0.0666 

6 

0.35875 

1.3937 

0.3937 

0.192 

0.07559 

7 

0.37125 

1.3468 

0.3468 

0.192 

0. 088156S 

8 

0.38375 

1.3029 

0.3029 

0.150 

0.045435 

9 

0.39625 

1.2618 

0.2613 

0.092 

0.0240856 

10 

0.40875 

1.2232 

0.2232 

0.044 

0.C098208 

11 

0.42125 

1.1869 

0.1869 

0.017 

0.0C31773 

12 

0.43375 

1.1527 

0.1527 

0.055 

0.0007635 

!  (w.  v  P  > 

_ : _ L_ 

=  0.3678949* 

_ 

! 

i 

* 

5 

*. 


Figure  1.  General  flow  diagram  of  the  aural  detection  model 
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Figure  A.  Horizontal  distance,  Ro>  from  source  to  shadow 

boundary  as  a  function  of  source  height  and  logarithmic 
sound  velocity  gradient,  3.  Height-of  observation  point 
is  5  feet. 


AVERAGE 


OCTAVK  RAND  LEVEL  Id*  RE.  0XW02  MtCROBAR) 


iwa>gj^s^eg>^&ieaajsaaq»<wegfejgBay^ 


~j-0-3  _j-°-3 _ 

ISO  ~30C  600  1200 

300  600  1200  2403 

OCTAVE  PASS  BANDS  (Hz) 


Figure  8.  5FL  to  senes  [4j 


.JIS 

S3* 

S3 

-JH 

-=3Sg 

Jll 


1 


3 

■33 


j 

.1 

-SS 


-  4$ 

m 


a 


- 

i 


Igg 

jft 


24 


2 

*i§ 

$ 


Preceding  page  blank 


Appendix  B.  COMPUTER  PROGRAM 

The  Gain  program  acts  as  a  calling  program  and  an  output  program, 
and  few  calculations  are  carried  on  in  this  part.  Most  of  the  critical 
task  calculations  are  carried  on  in  the  subroutines.  There  is  a  great 
deal  of  data  in  terms  of  arrays  representing  either  graphs  or  tables, 
and  a  large  percentage  of  the  data  is  read  in  by  assignment  statements. 
Besides  data  representing  graphs  and  tables^,  there  are  also  certain 
inputs  to  describe  the  conditions  for  detection: 

1)  DELTAT  -  value  of  the  temperature  gradient  (constant-5 F/1000 

feet) 

2)  TO  -  temperature  on  the  ground 

3)  T  -  temperature  at  height  of  target 

4)  DT  -  delta  time  -  time  incremental  value 

5)  WANO  -  wind  angle  measured  from  the  north  direction 

(y-positive) 

6)  FEET  -  height  wh^re  WIND  is  measured  (not  lower  .than  3  feet) 

7)  WIND  -  wind  magnitude  at  height  FEET 

8)  RH  -  relative  humidity 

9)  L  -  picks  out  terrain  type  (4  types  are  available) 

10)  XT  -  x-coordinute  of  target  position 

11)  YT  -  y-coordinate  of  target  position 

12)  ZT  -  z-coordinate  of  target  position 

13)  VX  -  x-coordinate  of  target  vector  velocity 

14)  VY  -  y-coordinate  of  target  vector  velocity 

15)  VZ  -  z-coordinate  of  target  vector  velocity 

16)  AX  -  x-component  of  target  vector  acceleration 

17)  AY  -  y-component  of  target  vector  acceleration 

18)  AZ  -  2-component  of  target  vector  acceleration. 

Discretion  should  be  used  when  inputting  target  positions  because 
there  is  a  time  limit  on  the  flight  path.  This  time  limit  could  be 
changed  but,  in  its  present  form,  120  seconds  is  the  time  limit  for  the 
aircraft  to  fly  into  detection  range.  This  is  true  uecause  the  aircraft 
does  not  necessarily  have  to  fly  into  detection  range.  If  it  dees  not, 
the  probability  of  aural  detection  will  not  reach  1.0,  and  the  output 
along  with  this  probability  will  be  "AIRCRAFT  DID  NOT  FIX  WITHIN'  DETECTION 
RANGE."  Therefore  this  time  limit,  which  is  changeable,  is  a  fact  that 
must  be  remembered  when  inputting  initial  conditions. 
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There  are  19  subprograms  in  the  model  plus  a  plotting  routine. 

Such  a  large  number  is  used  because  the  programmer  preferred  an 
uncluttered  main  program  and  an  orderly  assignment  of  tasks  to  a  sub¬ 
program  instead  of  grouping  several  tasks  to  one  subroutine.  The 
following  is  a  list  of  the  subroutines  used  and  an  explanation  of  their 
f unction: 

Subroutine  Called  From  Function  or  Contribution 

ACCl&iL 

MAIN 

Calculates  loudness  ( sones) 

ALPH 

Uf.T» 

Calculates  atmospheric  attenuation 

ALFKO 

SHADOW 

Calculates  residual  attenuation 

AL0CAT 

MAIN 

Allows  aspect  angles  (signatures) 

AVEAO 

ALPHO 

Locates  the  value  of  residual 
attenuation 

DEE BSE 

ACCUML 

Calculates  probability 

DIKT 

ALPK 

Mathematical  subroutine 

DVDTNX 

ALPH 

Mathematical  subroutine 

FLIGHT 

MAIN 

Calculates  new  target  positions 

FIG 

ALPHO 

Calculates  whether  an  observer  is 
in  a  sound  shadow 

INTEP 

FIG 

Interpolates 

INTER? 

AVEAO 

Interpolates 

ORDER 

ACCUML 

Ranks  orders 

RELPOS 

MAIN 

Calculates  elevation  angle,  azimuth 
angle,  and  range 

SHADOW 

"AIN 

Calculates  by  calling  other  sub¬ 
routines,  residual  attenuation 

TAB 

FUNCTION 

Allows  a  4-d.imensioned  array 

TABLE1 

ALPHO 

Table  look-up 

VEG 

MAIN 

Calculates  terrain  absorption 

WEED 

MAIN 

Calculates  phi,  the  wind-sound  angle, 
magnitude  of  wind  at  30  feet. 

Several  subrout 

'lies  were 

used  from  an  existing  model  [i].  The 

program  will  accept  target  signatures  for  five  frequency  bands  that  are 
critical  ones  for  an  airborne  source.  These  frequency  bands  are: 

75  to  150,  150  to  300,  300  to  600,  600  to  1200,  and  1200  to  2400  hertz. 
The  frequencies  used  in  calculations  are  the  geometric  mean  frequencies 
of  the  bands.  Velocity  and  acceleration  units  are  divided  into  vector 
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components.  Therefore  the  signs  of  the  velocity  and  acceleration 
components  must  be  remembered  in  relation  with  the  observer.  If  the 
target  is  at  positive  position  coordinates,  negative  velocity  coordinates 
must  be  inputted  for  the  aircraft  to  fly  over  the  observer.  The  sac?-' 
applies  to  acceleration  coordinates  also.  Figure  3-1  shows  how  the 
geometry  of  the  program  is  set  up. 


Figure  B-l-  Geometry  of  model  -  phi  is  the  wind-sound  angle  at 
height  30  feet.  Consider  the  wind  angle  wane  as  the  angle  of 
the  wind  vector  with  the  Y-axis.  Y-axis  corresponds  to  north, 
X-axis  east-west,  etc. 
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UQTDA.TaSA, 

ACCT(PN*D0N.P8C*3237O1O1OO»CC=77OO.OP=A2»JN*OOOO) 

FTN(L=0.PL=20000) 

LGO. 

EXIT. 

OOCODOOOOOOfcOOOOCeOOGOCO 

PROGRAM  MAIN  (INPUT. OUTPUT .TAP£5*INPUT.TAPE6=0UT?UT) 

DIMENSION  XtAR(50f  .YTAR{50)  .ZTAR'SC-' 

COMMON  ZAA{8)  *SIG{5»10.8)  « AH3 (5)  *D0 (10)  *C(6>  »EL  (6)-»YX(17)  »XK 
COMMON  FIGS (6.6) .2(6) .AC(5.4) .£(4) .EE(A» .SOMES (5.31 ) .K7ER.DELTAT,? 
COMMON  TG.PiE.rfAN0.XIND.FEET.PH.AX.AY.AZ.X7.YT.ZT.VX.YY.VZ 
COHKON/OIMCOM/TARDIM  H20J 
NftMEL I ST/OON/T ASDIM 

N  TELLS  HOW  MANY  SIGNATURES  ARE  TO  BE  READ  IN 
N  SHOULD  NOT  FXCCEO  JO-SEE  DIMENSION  STATEMENTS 
N=5 

TO  calculate  terrain  attenuation 
K  .  .  .  I-OPEN  area 

X  .  .  .  2-18  INCH  GRASS 

K  ,  *  .  3-LIGHTLY  TREED 

K  .  .  .  4-FOREST 

AT  IS  THE  TERRAIN  ATTENUATION  COEFFICIENT 
K=4 

KTER=K 

VARIOUS  CONDITIONS 
0FLT4T=-.5 

TFMRFRT»1R£  T  MEASURED  AT  2T 
T=38. 

TO=SO. 

PIE=3. 1415926536 
rfIND=50. 

FttT=15. 

VANO=?7C. 

RH=.40 
Ar=0. 

AY=0. 

AZ=0. 

VX=-338.98 

vr=0. 

VZ=0. 

KK  — PICKS  CUT  THE  CHOSEN  AIRCRAFT  FROM  SIGNATURE  SANK 
KK=1 

SIG(I.J.K)  I— PICKS  OUT  FREQUENCY 
J— A/C  TYPE 

K— PICKS  OUT  ASPECT  ANGLE 
READ  (5.21) 1((S1G(I rJ.Kl .N=l .5) » 1=1 .5) ♦ J-l .N) 

READ  (5.10)  (00(1) .I=1.N) 

READ  (S.21)  ( AHB ( I) » 1=1 .S) 

READ  (5.21)  (7AA(I)«I*l<.5) 

READ  (5.11)  (Cm.  1*1.6) 

READ  (5.11)  (EUI). 1=1.6) 

READ  (5.20)  <YX(I) .1=1.17) 

READ  (5.11)  ( (FIGS (I I »2J) »II=«»6) » JJ=1 .6) 

READ  (5.11)  (7(1), 1=1. 6) 

READ  (5.00N) 

READ  (5.21)  ((AC<I.J).I=1»5)»J=1.6) 

READ  (5s22)  (E(I). 1=1.4) 

READ  (5.22)  (FE(1).I=)»4) 

READ  (5.21) C (SONES(I.J) .1=1.5) .J=!»ll) 

10  FORMAT (RFlO.d) 

11  FORMAT (SFiO.O) 

12  FORMAT (7F10.0) 

20  FORMAT (8F10.0/8F10.0/F10.0) 


C 

c 

c 


21  FORHXT<5F10.0> 

22  FORMAT i*n  0.0 1 
DO  SO  X*1 *5 
ZTAR<kl*500«*K 
00  SO  J*1 «6 

YTAP  < J) *3280.84®  <J-1 1 

00  SO  1*1.8 

XTARd  1=3280.84®  <1-1  > 

XT=XTA»d> 

?T=Y7AP(J) 

ZT=ZTA3(K> 

CALL  WHOLE 
SO  CONTINUE 
STOP 
EM) 

SUBROUTINE  WHOLE 

OIMENSICN  F0<5) .A (5. 17.6) ,SPL<5f .SPLS15) .*0<5) .AT(S) ,AH<5*,AA{5) 
DIMENSION  OELADO(S) 

COMMON  ZAAf 81  «SIG(S, 10,81  ,AM9<5>  .COOO)  »C<6)  .Et<*J  .YXdTl  ,kX 
COMMON  FIGS <6.61 .Z<6) ,AC<5,4) »E <4) ,EEi4) »S0Nc3<5*il 1 »XT£R, DELTA! tT 
COMMON  T0,PI£«WANO.W!RD,FE£T.PH.AX.AY,AZ.XT,YT.2T.VX*VY»VZ 
COMMON/OIMCOM/TABOIH  <1201 
TIKE=C. 

DT*1 . 

L=0 

oo5i’i*i*s 

F=2.0*F 
I  FO ( I 1 =F 

OBSERVER  AND  TARGET  INFORMATION 
XO=0. 

YO=0» 

zO=5. 

CALCULATE  ATTFNUIATIOH  DUE  TO  HUMIDITY 

CALL  ALPH <SH. T. AH.FOl 
00  55  1=1.5 
DO  55  J=1 .17 
DO  S5  M=1.6 
55  A<I«J.K>=0. 

A<1, 1,11=1.5 
Atl,2, 11*3.3 
*<1.1.21*1.5 
A<1. 1.31=1.5 
AC1. 1.41*1.5 
A(l. 8,11*12.8 
Ad.*?, 11=13.3 
A(l. 10. 11=13. 6 
Ad. 11, 11=13. 8 
A<1. 12, 11  =  13.  *» 

Ad.  13,11=1*.  0 
A(l, 14.11=14.0 
A<1. 15,11=15.1 
Ad,  16. 11  =  15.5 

Ad.n,n=is.8 

Ad,ie.2)*A<1.10.31*Ad,10.41=li.4 
Ad.ll.2)=Ad.ll.31=AU.  11.  41=11.8 
A  d.  12.21  =A  d,  12, 3)  =A<i,  12.41=12.0 
AC} .13. ?1=A (1.13, 31 =Afl«13, 41=12.3 
A  <1. 14, 21  =A  (1. 14, 31  =Ad.  1^.41=12.4 
A <1. 15, 2) =A <1. 15.31 *A< 1.15.41=14.0 
Ail,16,21=Ad.l6.31*Ad  ,16.41  =  14.5 
A  d.  17.2)  =A  d.  17,3)  *A  <1.17,41  =  14.8 


i 
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A(2»l .11=4(2.1  *2)=.7 
4(2.2*15=2.5 

A<2.1*3)=A<2»1*4)=.7 
4(2*2, 31=4(2. 2»4}=.9 
4(2*1. 51*4(2. 1*6>=. 7 
A (2*5*6)=. 8 
4(2,10*1 1=16*2 
4(2«11, 11=15.0 
4(2.12*11=15.6 
A(?*J3«il=)5.8 
4(2,16.11=16.0 
4(2.15.1 1=17.0 
4(2.16*1 5=17.3 
A  {2*17*1 }=1 7 .5 
A{2.10»?)=5.4 
4(2*11*21=10,0 
4(2*12*21=10.2 
4(2. I3.?}=10.4 
A  (2« 14*21=10.5 
4(2* 15*21=12.9 
4  {2*16*21=14.0 
4(2.17.21=14.8 
4 52. 10*31=4 (2. 10 *4) =6. 2 
4  (2. H .31=4 (2,11 *4 )=3. 7 
4 (2* ! 2.31 =4(2* 12*4 1=7.0 
4(2*! s.3l=4(2»13*4}=7.4 
4(2. 14.31=4 12. 14*4>=7.S 
4(2* 15.31=4 (2.!5>4>=10. • 

4 (2* 16*31 =4 (2. 16*4)=12.0 

A  (2*17.21=4  (2.  17«>.1=13.0 

*(2.10.51=4.8 

4 (2* II *5! =5.3 

4(2.13.51=6.0 

4(2,14,51=6.1 

4(2.15.51=9.8 

4(2.16*51 =10.3 

4(2*17*5) =11.8 

A  f2* !?. 51=5.8 

4(2.10,61=4.0 

4(2,11 .61 =4.4 

A(2. 12.61=5.0 

4(2,13*61=5.2 

4(2*14*61=5.3 

4 (2* . 5.61 =8*0 

4(2*16, 61=9. 8 

4(2*17, 6)=>1.0 

4(3*10, 11=16. 6 

4(3*11*11=17.? 

4(3*12*11=18.0 
6(3*13.11=16.2 
4(3.14.11=15-4 
4(3.15,11=19.3 
4(3,16,11 =20 . 0 
A(3*17. 11=21.0 

-*(3*10T£'!=4(3.;O,3)=A(3.1O»4  5=fi.O 
4(3.11 .21 =4 (3* 11, 3} =4(3*11 ,41=8.6 
4 (3, 12.?) =4 (3, 12.3) =4 (3. 12*6) =9.0 
4(3*13*?1 =4 (3,13*31 =4(3,13*41=9. 6 
4 (3. 14,?} =4(3,16, 3) =4 (3, 16. 4) =10.0 
4 (3,15,21=4(3.15, 31=4 (3* 18, 4} =12.0 
A (3* 16.?) =4 (3. 16, 3) =4 (3, 16, 45=16.0 
A'3,1 (,?;=A(3.17.3)=4{3»i7,4)=14.P 


A ( 1 *5* 1 ) *6. 1 
A(l«6»l)*ll*3 
A(l,7tl)*12.0 
A ( l *2*2) *2.5 
A ( 1 *3*2) *3,0 
A(l*4*2)*3.8 
A ( 1 *5*2) *4.5 
A ( 1 *6*2) *9.0 
A (1*7*2) *9.9 
A(l*8.2)*10.25 
A ( l *9*2) "10.8 
A ( 1 *2*3) *2.5 
A ( 1 *3*3) *3.0 
A ( 1 *4*3) *3,8 
A ( l *5*3) *4,5 
A ( 1 *6*3) *9,0 
A ( 1 *7*3) *9,9 
A ( 1 *8*3) *10,25 
A ( 1 *9*3) *1 0,8 
A ( 1 *5*4) *3,9 
A ( 1 *6*4) =4,5 
A ( 1 *7  *4) =5,0 
A ( 1 *8*4) *5,5 
A ( 1 *9*4) *6,0 
A(l*10»4)=6,25 
A ( l ♦ 1 *5) *1  *5 
A ( 1 *2*5) *1 .6 
A ( 1 *3*5)*1 ,7 
A ( 1 *4*5) *1 ,7 
A ( 1 *5*5) =1 *8 
A ( 1 *6*5) *2,0 
A ( 1 *7.5) *2,6 
A ( 1 *8*5) *3,2 
A  ( 1 *9*5) *3,8 
A (2*3*1 ) *4,0 
A (2*4* 1 ) *4,8 
A (2*5*1 ) =5,9 
A (2*6* 1 ) =9,2 
A(2*7*l)*10.6 
A(2*8*l)*12.0 
A  (2*9*1 ) *13,2 
A (2,3,4)*1 ,0 
A (2»4*4)=1 ,5 
A (2*5*4) *1 ,7 
A(2.6*4)=2»5 
A (2*7  *4) *3,2 
A(2»6*4)*3,9 
A (2*9*4) =4,2 
A (2*3*2) =2,0 
A (2*4*2) =2,35 
A (2*5*2) *3.0 
A (2*6*2) *6,9 
A (2«7«2)*7,7 
A (2*8*2) *8,25 
A (2*9*2) *8.9 
A (2* 1 *5) **8 
A (2.2*5) **9 
A (2*3*5) *1 .0 
A (2*4*5) *1 *0 
A (2*5*5) *1 ,0 
A(2»6*5)*l *7 
A (2*7*5) *2,4 
A (2*8*5) *3.0 


A (2*9*5) *3.65 
A(2«3*3)*1.0 
A(2»4«3>*1 .5 
A(2*5*3)=2,0 
A(2*6*3) *3.8 
A ( 2*7*3) “4.5 
A (2*8»3)=5.2 
A(2.9»3)=5,8 
A(3*l*l)".8 
A(3.2*l)=2,2 
A(3«3«l)*3,3 
A(3*4«l)=4.0 
A (3*5 . 1 ) =4,6 
A ( 3*6* 1 ) “10,5 
A (3*7* 1 ) *13,4 
A(3*B.1)*15,0 
A(3*9»l)*16.0 
A (3.6.4)*3,1 
A (3*7*4) =4,0 
A (3*8*4) =4.8 
A (3*9,4) =5,4 
A (3*1 *5) *.8 
A (3»2*5)=*8 
A (3»3»5)*«9 
A (3*4*5) =1 .0 
A (3*5*5) =1 ,0 
A (3*6*5) =2,0 
A(3.7*5)=3.0 
A (3*8 .5) =3,8 
A (3»9.5)=4,3 
A (3, 1 *2) =*8 
A (3*2*2) =1 *0 
A  (3»3*2)=2.0 
A (3*4*2) *2.4 
A(3*5*2)=2.6 
A (3*6*2) =4 ,8 
A(3*7*2)=6.0 
A(3«8*2)=6,8 
A (3.9*2) =7 .4 
A(3.1*4)=,8 
A (3,2*4)=1 .0 
A(3*3.4)=l,2 
A(3»4*4)=1.8 
A (3.5.4) *2,0 
A (4*5* 1 ) =4,0 
A (4*6  *  1 ) *9,0 
A(4*7*l)*10. 
A(4«fl*l)=ll»0 
A  (4»9*1)=11  ,6 
A(4*5*4)*3.6 
A (4*6*4) =5 . 0 
A(4*7»4) =5.9 
A (4*8*4) =6,5 
A (4*9*4) *7.3 
A (4»5»2)*4,0 
A(4.6*2)*9.0 
A (4. 7* 2) =10,0 
A (4*8*2)*1 1.0 
A(4*9*2)=11.6 
A (4*5*5) =1 ,8 
A (4*6*5) =2*8 
A (4*7*5) *3,9 
A (4*8*5) =4,6 
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,  ww*w.»->!w.v?v  (■**** 


c 

c 


c 

c 


A{4.9*5)=5.4 
A,'4.5«3)=4.0 
A (4*6*31=9.0 
A  (4*7 *21=10*0 
4(4,8.31=11*0 
A {4.9*31=1 i .6 

4<5.1.!)=A(S.1,2)kA(5,1.3)=2*C 

A(5.2.1)=A<5.242)=A<5,2r3>=3*0 

A <5.3. 1 >=A <5*3. 2)=A (5*3. 31=4.0 

A (5 *4*1) =A <5*4, 2) =A (5*4*31=4.6 

A(S*5*1)=A{5»5*2)=A(S*5*3)=5.5 

A{5.6*1 1=2 (5*6*2) =A  (5*6*3)  =12*4 

A (5, 7*1) =A <5*7*21 =A <5 *7.31=14. 8 

A  (5. 8* 1>=A <5*8*21=4 (5*8. 3) =16.0 

A<5.9.t)=A{5.9*2)*A(5.5*3)*l6.4 

A(5tl0«i')  =  A{5«10.2)=A<5*10*3)  =  l8.0 

A(5,S*4)=5.8 

A{5»6*4)=7.0 

A (5*7*4! =8*2 

A (5*8*4) =9*3 

A<5*9.4)=10*0 

A<5. 10*41=10. 8 

<t(5, 1*51=2.0 

A(5*2*5)=2.2 

A(5.3*5)=2.4 

A (5»4*5)=2»6 

A(5*S. 51=2.8 

A <5*6»5)=4,0 

A (5*7*5) =5.2 

A(5*8.5>=6.0 

A<5*9*5}=6.8 

t <5*10*51 =7.4 

KKK=0 

H=0 

STARTING  POINT 
VRI7F  (6.802)  XT.YT.ZT 

2  TTME=TIME*DT 

CALL  WEED§WANO,tfTND,FEET.XT»YT*ZT*XO*YO.ZOrPHl,W30> 

CALL  PELPOS{XO*rO,ZO.XT.rT.ZT«ZA.flETA.O.VX*vr,VZ> 
IF(RETAcLT.2.)GOTO  45 

#«****»*«*# 

CALCUALTE  attenuation  DUE  TO  terrain  type 
K*KTER 

CAls.  VEG(BETA.O*AT*KvAC*E»EE) 

*******«********>*#*****-»**«*e*  .•e&*3»*«**4****«*#«**»*».9*****o«a**»o 

XX=X7-X0 

CALCULA  E  ATTENUATION  DUE  TO  TU98ELENCE 

CALL  5H(  J0W(XA.A«RETA.EL»D.DELTAT.TtT0*FIG5*A0'PHTc'*(30»XT.YT.ZT.S0 
l.YO.ZO.YX.C) 

IF  <  AO  (ll.F.O*  <-!.))  GOTO  45 
00  3  1=1*5 

3  AA(I)=AH(I)»(D-D0(KKJ  1/1000* .AOn-.ATd) 

00  90  1=1*5 

CALL  ALOCAT  <  I  *>(K*SIG»ZA»ZAA  .REWARD) 

SPLS (I) =REWAR0-2f .*AL0Gjfl (D/DO (KK) 1 *3.-AA ! I ) 

90  C0M7INUE 

9*4*m*t*****««*«(4»**««**«*«««**t««:*«t«4»«)*«ii**ii««****i)««t)jiit 

COMBINE  SPL  WITH  AMBIENT 

4(l«*4C4**44<4**e*4<!'444r4*H4«****(et4»*»*e4*4**«»*«4*4«4<  ?<•£«••  9 

00  1000  1=1 45 

SPL<i>=io.**(sPLsm/io,)*io.«*<AMB(r)/io.) 

SPL ( 1 1 =10  **ALOGIO (SPL (11) 
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OGO  CONTINUE 

CALCULATE  PROBABILITY 

CALt.  ACCUHL  (SPLi  AH5*S0Nf  S«VAl  UF  .ARC! 

IF(ABC.ES.C.)50Tfi  ^5 
L*L*  I 

35  IF'L/61*L.F0.L>WSITE  <6.8005 
>F'L/6l*L.F0.L)L-0 
Wo-TE  (6*23^  TIME.  VALUE 
IF  CVALVt  *SE  .  ,9R99<?9999>  60T0  19 
45  CALL  FLTGHTU:<.AV,£T,yx,yy,VZ,XT«YT.ZT,DTl 
IF  <  T ! HE .GE  » 1 20  » ) GOTO  SO 
GOTO  ? 

15  CALL  tfELP05(XC»,YG,ZO«XTtrT*ZT,ZA.8ETA-O*Vi;,VY.VZ; 

KKK=2 

WRITE  (6.28) 

WRITE  <6»?3)  XT.YT.ZT 
XPITF  (6.74)  BETA 
WRITE  (6.400)  ZA 
WRITE  ( 6  ♦  7  0  ’•  U30.PHI.TIME 
WRITE  16.100) 

KsKTER 

IFfX.FO.DVPITE  <6.755 

IFi':c.SG.2iWRITF 

Ir  iK.E0.5)WR!TE  56. 7?) 

IF  CK.F0.4')  WRITE  (6*73) 

WRITE  <6*)4)  {FO(j ) ,  AH  <  I ) *AC ( I) *AT (I) »I=i .5) 

WRITE  «6.1Cf)) 

WRITE  (6.16)  (EOSJI.AAI J)«SPL?.iJ)«AM8(J).J=l,S) 

WRITE  (6.29)  ft 
WRITE  (fc.100: 

XD*X7 

Vr.-yr 

ZO=ZT 

VS*1080.*SO»T( (T*acr»»/(T0*459,>  > 

OEf.AY=r./VS 

CALL  FLIGHT (Ax . AY, £Z »VX . VY.VZ.XD. Y0.2D.0EL AY) 

CALL  RELPOS « AO  .  Yft  *  ZO  .  XD  »  YO  .  ZD*  Z  A  *  R£  T  A  ,  f>0  ,  VX  «  VY  »  VZ  > 

WRITE  ‘G.ftOD  DEL  Ay, 00 
DT“9  - 

TJHFrTtHE.DT 

c-mjl  flight  (Ax.  ay,  az«vx.vy,vz«.xt,yt»zt.dt) 

CALL  RFLROS ( 20. YO. 70. XT .YT.ZT.ZA.SF,  T A  »D»VX»  YY * v,*2 ) 

WRITE  (6.100) 

WRITE  (6.73)  XT.yt.ZT 
WRITE  (6*74)  R£TA 
WRITE  -:6,4-jO)  ZA 
WRITE  (6.15* 

SO  COM TJIJIjE 

IFJKKK.EO.OVRtTF  (6,200) 

14  F0r'HfT(lX,6H  FREO=,F10.R.2X,aH  AH=.Ffs.3.;12w  OR/, 000.  FT.2X.4H  AOe , 
6rAs3,?Xr6«  ATr.FR,3) 

15  rna MaT.-lA.3SK  VISUAL  DETECTION  OF  TARGET  AT  TIK£,C8.2>9H  AT  RANGE, 
6F13.2) 

16  FGPPArdX«6H  tR3?=.F10.?,?X,i5.-i  TOTAL  ATTEMUATI 0N=  .FR ,3 ,2k,SH  S?L  = 
6  »  F 1 0 . 3  »  2  *  *vH  am3 I fnt s»F,0»3) 

?3  FORM  A  T (1 X , f  7 . ? , 1 1  X , F 1 0 . 5 ) 

?8  FORMAT  ;  mi ) 

?4  FORMAT  (1X,31H  OBSERVER  HEAPS  TARGET  t.  7  RANGE, HO.?) 

70  FORHATUX,SH  U30=,R6.2,C-H  C<hj..,f6.2,6W  T lMF=tF#,s;>*, 

73  FORMAT ( I X.4h  vT=.FP.?,-H  YTt.fp.^.ah  ZT=. F8. 25 
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74  roRMAT<]X.6H  5ETA=«F6.2) 

75  F ORH  AT ( 1 X « 1 fiH  T£P»A5N-OPEN  AREA) 

76  FORMAT UX«?2K  TERRA IN- 15  INCH  GRASS) 

77  FORMAT (1X»22H  TERRAIN-L  IGHTLV  TR£FO> 

7g  format hx»15m  terrain-forest} 

100  FOFMAT<}X«//F) 

200  FORMAT { 1X»42H  AjPCRAFT  0!t>  NOT  FLY  WITHIN  HFARI NO  RANGE! 

400  FORMAT (1X»1SH  AZIMUTH  ANGLE=*FI0»25 

600  FORMAT { 1M1  *"X-*5H  TIHE»5X»'54K  PROBABILITY  OF  AUDITORY  DETECTION? 

801  FORMAT { 1X« 14H  DELAY  TIM-  IS »F4,6«32H  GrCONOSt  TAPGfT  IS  NOW  AT  RAN 
SGF.F5.2) 

602  FORMAT UMi.2X>5H  T*ME«SXt.34H  PROBABILITY  OF  AUDITORY  DETECTION* 1  OX 
6*1«H  STARTING  P0?NT*3X«4H  X.7=*F8*?«4K  YT=*Fft.?e4H  ZT=,F8,2) 


RETURN 

END 

SUBROUTINE  ACCU«L ! S?L  *  A MB »  SOMES ♦ PP08  *  ARC 1 

DIMENSION  50NFS!5*11! * SPL 15) . AK£J5} »S<S> *SS(5) .OECfll 5 

DATA  0EC/20.«10*»40*.»S5.«S0..75.*8e.»Re.»lO‘>*.  HO. *120./ 

ASC=1 * 

SAH3=SNL=0* 

00  25  1=1*5 
K=6 

:f«SPI.( I) .LT,2«*.OR.SPL(I).OT.120,)GOTO  20 
no  25  ..<=1*10 

I F  { SPL  n  > . GE  .DEC  f  J)  *AND. SPL  <  I )  .L? ,OEC  I J*  I  > )  GOTO  1 0 
GOTO  15 

10  Sm  =  <(SPLU)-OEC(.J))/(OEClJ»)>-DECUm*fSONESCl.J*I>-SONES{J»J>) 

5*S0NESa«U> 

I "  continue 

CALL  ORDER  IS! 

GOTO  »5 
20  A?C=0. 

RETURN  -~  .  '  - 

?5  CONTINUE 
00  35  1  =  1-,  5 

IF ! AMR ( J  J  »Ll .20.. 09.  AMS  <n .G7.!?0.)GOTO  20 
DO  35  J=1 *10 

IF (AMR. I>  «Gt • OEC(J) .AND* AMR  II) »LT«.OECf  J*1 J I  GOTO  Vi 
GOTO  3= 

30  S5 { I )  =  i  { AH??  { I )  -DEC  <  J)  5  /  fOEC  ?  J*  1 !  -DEC  ( J)  J  ?  ♦* ' 50N€<-»  I  t  J*  1  *  -SOM*  S ;  I  *  Ji 
6i*S0NES«I«J) 

35  CONTINUE 

CALL  GPDERSS) 

00  40  1=1*5 
S4MB=SS-:I  s«SAMB 
SNL=S 2 1 ) *3NL 
40  CONTINUE 

CALL  0EE8EE ?SnL»SAMB*R90Bt ABC) 

RETURN 

END 

SUBROUTINE  ALPH(Rh,T,Am,FO> 

DIMENSION  Pf  17»11)  «W(17,n>  ,FOlS)  «AHI5)  .TTU7!  *F!  <5) 

55=20.3 
00  i  1=1*17 
IT  <I)=8B 
1  8R=B8»S.O 
00  2  1*1 *17 
P;j>i>*,l 
P'I»2)».15 
Rc!»3)=*2 
P{I*4)=,3 
Pi  1  *5>  *■»* 
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R(I«6>=.5 
P(I.7>*,6 
P(I»S)*,7 
R(I*9)n.8 
9(1.10)*. 9 
2  R (I . 1 1 )=> ,0 
an*i;=.3 
W‘1.2)*-,-42 
Iff!  *3)=.57 
wf»«4)=.es 
W(i*S>=i.!2S 
¥(! 

W { 1 «?}=1.7 

wn,s)=i.9 

W ( I  ♦  <5}=2»125 

W(l,10i=2<*5 

W<1»H)=2.7 

V (2. 15 =.36 

V(?*2)=.53 

V {?  >3)=.7 

W<?,4)=1.1 

«(2*5)=1.45 

W(2,6)=i.70 

V(2^?)=2.1 

V(2«8)=2.46 

W(2.9>=2.80 

W(?«10!=:.l 

Wf2.ll ,=3.5 

W(3.1)=.4fc 

W(3»2)=.68 

W(3»3'=.*° 

Wf3.4:-=1.3 

Wf3»5»=l .75 

W(3»6>=2.1 

Wt3.7)=2-fe 

W*3.5)=3.G5 

W{3»9)-3»5 

W ,3.10)=3.9 

W(3.!i;--4.3 

W(4,})=.55 

W(4,2)=.«l 

W  (4.3)=1 . 1 

y{4.&>=l.6 

V(4,5)=?.07 

V  (4.6) =2. 7 

W(4.7}=3.3 

Wff'*R)=3.8 

W<4»9)=4.2 

V(4,]0)=4.6 

W<4,11)=5.2 

W (5, 1 ) =.f 7 

WfS*'Ms,9R 

W(R.35=1,3 

<j(5»4)=1  .9 

W.'S,5)=2.6 

vi  (5»8>  =3.2 

WfS.Tjs:  1.9 

W (5.8) =4.5 

W (5*9) =5. ! 

W(-5.i3)=5.9 

W*5. 1 ! )=6.» 

Wf6ii >=.B 

W(6.2J=t-2 


W(6»3>=1.5i 
V(6*4)=2.4 
W(6.5)=3.2 
W(6.fe)=3.9 
W(6<7)=4.8 
V(6-*8)=5.5 
W(6«9)=fc.2 
W1:*.10)=7.p 
f  W(6»ll)f*7«8 
W<7.1)--.94 
'J  t  7c2)  >*1  .4 
Wt?.3)=i .8 
WC7.4J  =2.8 
W<7,5>=3*7 
W(7 .6) =4.7 
W{7 .7) =5.6 
W(7.8:=0.5 
W  (7-.9)  =7.S 
*(7.10)=6.4 

>»  ( 7.1 !  )  *9.} 
W(8«l2=1.2 
W(8,2)«l .65 
W <8* 3) =2.2 
V(8.4)=3.4 
W(8,S)*4.5 
W(8.6)=5.5 
W(S.7)=6.6 
W(5.8',=7.R 
W15.9>=9.0 
W(8«10)=10.0 
W(8. 115=11.0 
W£4»l)=1.3 
WC9*?»=l.o 
W(9,3)=2.6 
W (9.4) =3.9 
*(9,55. -5. 2 
W(9.6)C6.5 
W(9.7>=7.0 
W  (9»R)=9. 1 
V(9.9)rlO,5 
W(9il(5j=12.0 
W(9.ll)=13.0 
W(10. 1!=1.6 
wn%?>=?.4 
W(10.3)=3.1 
W( j0.4)=*t8 
WU9.5>=6.2 
W(!f).8/=7.9 
W ( 1 0. 7>=9»2 
I  V(10. &)=!!. 0 
WU0,9)  =  12.5 
Wfl0,10)=l4.0 
w < 1 0.1 2  >  =i5.5 
W(11.1J=1 ,8 
W( 1 ! .2) =2. 7 
W ( 1 1 .3) =3.6 
Wt!J »4J=5.5 
=fr.!'5)=?,3 

w<n.6?=^*o 

W(1 J .71=) 1.2 
W(ll .81=12.7 
W(H»9j=l4.5 
W(1J. 101=26.0 


W « 1 1  *  111=18.0 
W(12.1)=?.3 
W (12.2) ~3»3 
W ( 12.31=4.3 
W* 22,4) =6 .5 
W (12.5>=A*7 
W(12»6>=12.2 
i  W(l?. 71=23.0 
1  W (22*8) -25. 9 
Wa2.9)=!>.> 
Wd2,10»=f9«0 
1  W(12.2i)-22.Q 

Wcl3.1)=2.5S 
W<)3.2)=3.8 
W(23.3)=5*0 
W(13t4)=7.S 
W(23.5:=20<0 
W ( 23.65=52.5 
W(13.7)=;5.() 

W ().3»8)  =27.5 
W (13.9) =2<5«0 
W(13«10)=?2«5 
W(12fH)=2S.O 
W ( 14*1 ) =3.0 
Wi)4,2)=-l.5 
w;i4.3)=j-9 
W ( 1 4  *4) =8.9 
W ( 14.5) =12.0 
V<14,6)=!4.5 
W(14<7J=17.5 
W ( !4c8) =20.5 
V(14.9)=23.5 
W(14. 10=26.5 
W(14.11 )=29.5 
W (15* 1 )=3.4 
W<2S,2)=S.O 
W(1S.3)=6.7 
V ( 25,4) =10.0 
W»13t5)=!3.5 
W(!5»6)=17.0 
W il5»7)=?0«0 
W ( 15.8) =23.9 
W(15.9)=27.0 
V(15*2P)=30.0 
W (l5rl 1 1=33.0 
W(I6»1)=4.0 
W ( 18.2?  =5.9 
W  f 18,3) =7.9 
SVlo,4>=12.C 
W(16«5)=16.C 
'Jj,l6.6)  =  19.5 
W(16.7)=24.0 
W(18,8>=-28.0 
W(16,9>=3?.0 
W(18»10)=35.0 
W  (16.! } )=38.0 
ic{17rH=4.S 
W{]7.2)=£,8 
Wd7,3>=9.0 
W { } 7»4)=23,5 
17 .51=18.0 
Wfl  :.6)=23.C- 
I  W  f 1 7*7) =27. S 
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W(17«3)=32.0 

¥(17.93=36.0 

¥{17,103=M.0 

¥(17,113=45.0 

CALL  DINT  (R.W.TT.l/.ll-T.RH.WW) 

DO  50  K=l**» 

FI  (K>=FC(K>/1G00. 

AH(K}=(,1«F1  (K)*(T*45.) )/( <F1  {K3/WW**23*((WWA*2)/f1  (K3>) 

RO  continue 

RETURN 

END 

SUBROUTINE  AlPH0IXX,A.6Fr4.FL»D»0ELT4T,T.T0,Z0»ZT.FIG5,A0,yX.Z.PHI 

6»3*c*voo.xt,yt«xo,yo3 

DINEN^rON  Z (6) .C<fc) »FI65<6.63 »EL (63  «YX(17) .4(5,17.63 
DIMENSION  GIG (5» 10.3} «AM3<53 ,AQ(5> .0£LAD0(5>  .AH(S) .00(10)^(53 
COM«ON/OIMCO«/TA3DIN  {120 

so  format (57H  observer  is  in  the  shadow  region  at  tkf  given  time  belc 

6W* 

R=S0RT  C (XT-X03  **?♦ { YT-YO) **2) 

5  CONTINUE 

J.l=l 

IF(RCTA.LT.,2«)G0T0  3G 
IF{BETA.F0.2»)GOTO  10 
IF (BETA.GT.2. . ANO.BETA.LT* 10. 3  GOTO  15 
GOTO  20 

10  00  11  1=1  ,5 
NN=I 

call  a ve ao  t  xx , a , be  ta . fl . q .nn , aqb ar . yx  > 

CALL  TABLE! (i.O.SETA.BBT.BR*) 

K  (BBT .EQ*1 «0.0R.10BAR»F0»0»  3G0T0  20 
AO(l3=&03AD-.3*B87»{DELTAT.4.S)*i12*8B!<*(W30-R.,>) 

11  CONTINUE 
35  NM=! 

CALL  4VEA0(XX.A«8ETA.fL»D.NN.A0BAR.YX) 

CALL  TABLEIU.O.BETA.BPT.BBW) 

IF(BBT.FO.l.G.GR.A0BAP,F0,C.iG0TG  20 
40n)=A0RAR-.l*BRT*{0FlTATo4.5}+.12«»9BW«fW30-c.S> 

JJ=? 

20  00  ?c  J=JJ*5 
*iK=J 

CALL  AVEA0(XX,A»PETfi.F{.,0,N4«*,A0BAR»YX3 
:F<A05*o.E0.0.3GOT0  30 
A0(.:)  =AOBAR 
??  CONTINUE 
RETURN 
ENTRY  AL 
.3.3=1 

R=SC1RT( (XT-XO) **?♦ (YT-YO) e»2) 

CALL  FIG,B«C»ZT*Z*P0«“IG5) 

IF (RO.EO. 0.3 GOTO  5 
IFtR.Lt .R03GCTQ  5 
IF(S.G£,3.*»OJGOTO  25 
GOTO  27 

25  IF ( PHI. LT. 110. -OR* PH! .GT. 1 50. 3 GOTO  27 
WRI7F  (6-S03 
00  26  1  =  1,5 
40(13=30. 

?fc  CONTINUE 
?9  CONTINUE 
SFTURN 

27  00  23  11=1.5 
NN=!I 

CALL  AVEAO (XX » A.BETAjFL .RO,N?«,  408AR,  YX3 


I 


4 
* 
'  I 

-TS 


1 

3 


-1 

# 
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IF  (AORAP.EQ.O. ) GOTO  30 
A0(II)=AQ8AR 

OELADO(fI)  =  (R-RO)*  (30.-40'.  II) )/(2,»R0) 

28  AO(II)=AO(II)*OELADO(!n 
WRITE  (6*50* 

RETURN 
30  AO(l)=-l. 

IF (D.LT >6000 • /GOTO  35 
RETURN 

35  no  40  1=1*5 
40  AO(I)=0. 

RETURN 

END 

SUBROUT I N£  ALOC AT ( I . KK ♦ S I G . ZA . Z A  A . RE WARD ) 

DIMENSION  SIG<5.10, 8), ZAA(8> 

IF(74.GT.180.)ZA=360.-Z4 
DO  5  N=l,4 

IF (ZA.GE.ZAA (N) . AND. ZA.LT. 7 AA (N*l  > )KKK=N 
5  CONTINUE 

IF  17A.EO.180. )KKK=4 

f : uZ2,At  (KKK) > ' {?AA  >  -zaa  (KKK  j  1)  *  (Sir,  ( i  .kk.kkk.  i  »  -slot 

j I . KK .  <K|< (  J  *SJG( I. KK.KKK) 

RETURN 

END 

SUBROUTINE  AVEA0(XX.A,BEIA,EL.D.NN,A03AR,vX) 

DIMENSION  EL(6)«YX(17)»A{5,I7*6) 

IF(D.lT.GOOC.)GOTO  20 
IF (BETA.LT. EL (1 ) >GOTO  20 
DO  10  N=l,5 
I=N 

IF (BETA. 6£. EL (N) .AND. BETA. LT.EL(N*J))GOTO  30 
7  CONTINUE 

If (BETA. EO. EL (6) )GOTO  30 
20  CONTINUE 
AORARsO . 

RETURN 

30  0K=0/1 000 . 

DO  40  K=i.JS 
J=K 

IF (OK.GE.YXiK) .ANO.DK.LT.YX (K*l ) )GOTO  50 
40  CONTINUE 

IF{nx.FQ.YX(lG))GOTO  50 
GOTO  20 

50  CALL  INTFRP(OK,I.J.A,nN.BFTA,EL.Vai  UE.YX) 

AOBAR=VALUE 

RETURN 

END 

SUBROUTINE  DEFPEE ( SNL . S AMB . PROB .ARC) 

DIMENSION  SL (]R; .OJi (JR) ,012(10} 

DIMENSION  P ( 1 2) .BAND (13) »BANHID(I2>  »P0(12) 

)0A«5  SOrt10o^!5|*«ft°y’?5*',30”35**40**A5**S0*’SS*'60*,65.»70..75,.80 
J • »G5. . 90. *95. .100./ 

DATA  Of 1/.76,. 47.. 34,. ?o,. ?7,. 27.. 27.. ?7.. 27.. 27.27., ?7». 27,. ?6,. 2 

DATA  Q^/  et,  52.  ^,.32°. 32^32,  >324  >32^3?^3?^32^  32<  32 

jJcwJf  ♦.Jl  * #27/ 

^DATA  °/» 005.. 91 7.. 044.. 09?*. 1 50*. 1??.. 192.. lie. .0R2,. 044.. 01 7,. 005 
RRGB=0. 

Ir (SAH8.LE,0. .OR.SNL ,LF.O.)GOTO  25 
T4M9= (6LOG10 (SAMP) .1 .?> /0.03 
SSS= ( AL0G1 0  <SNL) *1.2) /0. 03 
IF (TAKB.LT. 15.. OR.SSS.LT. 15.) GOTO  25 
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DO  1  1*1. IP 

IF (TAMB*GE»5L { I ) • AND, T^MB.LT.SL (1*1)) K* I 

I  CONTINUE 

IF (K,FQ»0)GQTO  5 
XBAR1= (QI 1 (K) *012 (K) )/2» 

XRAR2=(0I1  (K*l)*QI?(K*l))/2„ 

XBAR=  (  (TAMg-SL  (K>  > / CSL  <K* 1 )  -St.  (K )  )  )  *  { XBARP-X8AR1 )  .XBAR1 
A= { (TAHB-SL (K) ) /(SL (K*l J-SL  fK)))*(QIl(K*l) -Oil (K) ) »QI 1 (Kl 
B=( (TAHB-SL(X) >/(SL<K*l)-SL(K)))*(QI2(K*))-0I2(K>).QI2(K) 
STC,HA=SORT  { <  A-X9AR J  **2*  (9-XBAR) **?) 

BANfHI)=X8AR-3.«SIGHA 
00  2  1=2.13 

?  BAND (I )=BANO < I — 1 > *SIGMA/2« 

00  3  1*1,12 

3  9ANMID(I»  =  <8AN0n).BAND{I*l)>  ' 2 . 

0ELTA=SSS-TAM8 

DO  A  1=1,12 

PD ( I ) = fOELTA/BANMIO ( I >  > -1 . 

IF {PO (I) ,LT. 0 , ) PO  ( I ) =0 ,0 
IF(PD(I).GT.1.0)PD(I)=1.0 

4  PR0B=P0(T)*P<I)*PR08 
IF  (PROS.LE..O.  )  PROB=C. 

RETURN 

25  PROB=0. 

RETURN 

5  ARC=0. 

RETURN 

ENP 

SUBROUTINE  PINT  (R.W.TT. I . J.T.RH.WW) 

DIMENSION  R(I.J).  W(I,J).  TT(I) »  X(100).  r(100) 

P(I.J)  CONTAINS  INDEPENDENT  VARIABLES.  VALUES  MUST 

HAVF  SMALLEST  value  AT  J=1 

Vtl.J)  CONTAINS  DEPENDENT  VARIA8LES 

TTU)  IS  THE  PARAMETE®  THAT  DETERMINES  WHICH  OF  THE 

I  TH  CURVE  OF  P(I.J)  AND  W(I.J>  TO  C0N5IDFR 

T  15  THF  KNOWN  PARAMETER  VALUE 

RH  IS  THE  KNOWN  INDEPENDFNT  VARIABLE 

WW  IS  THE  OEPFNDFNT  VALUE  TO  BE  SOLVED 

IF (TT (1 ) -T) 1.3.2 

1  IF (T-TT {1)19.6.5 

2  WRITE  (6.17) 

3  DO  4  TI=1,J 

x<in=Pti.n> 

4  Y(IIt=«(MI) 

goto  a 

5  WRITE  (6. ! R) 

6  00  7  11=1. J 

xni)=F(i.n> 

7  Till) =W (I. II) 

a  CALL  DVD I NT  (AH.WW.X.Y. J.2> 

RFTURN 
R  L=2 

10  IF'TT (L)-T)11«12.14 

II  L'-L  ♦  1 
GOTO  10 

12  DO  13  11=1. J 

xnn=R<L»iJ> 

13  T(in=K<L.II) 

GOTO  8 

14  LL=L-1 


00  15  11=1. J 
X ( 1 1 )=P {LL . 1 1 ) 

15  Y(II)=W(LL.II) 

CALL  0V0INT  (RH.Z.X.Y. J.2) 

00  16  11=1  * J 
X ( T I ) =P (L» I T ) 

16  Ymi=WIL.II> 

CALL  DVDINT  (RH.ZZ.X.Y.  ).?> 
WW=Z*{ZZ-Z>*(T-TT<LL) )/<TT{L>-fT(LL) ) 
PFTURN 

17  FORMAT (20H  VALUE  OF  T  TO  SHALL) 

IS  FORMAT (20H  VALUE  OF  T  TO  LARGE) 

EK'D 

SUBROUTINE  OVDINT  (X.FX.XT.FT.NP.ND) 
DIMENSION  XT(NP).  FT (NR).  TT(16) 

N=ND 

Nl=(N-l)/2 
N?=N/2 
N3=NP-N2*1 
IF (NP-N) 18*1.1 

1  N4=Ni*2 

Ir (XT (1 ) -XT (2) >2.26*19 

2  CONTINUE 

IF <X-2.«XT(1).XT(2> ) 15.15.3 

3  IF  (X-2.*XT (NP) *XT (NP-1 1)4.4.15 

4  IF(NP.LT.10)G0TC  6 
N5=NP-N 

5  N5=N5/2 
N6=N4*N5 

IF (XT (N6) .LT.X)N4=N6 
IF (N5.GT . 1 ) GOTO  5 

6  TF(X-XT(N4) )9,7.7 

7  jF(N4-N3)8.9.S 

8  N4=M4.) 

GOTO  6 

<S  N4=N4-1 
N5=»(4-N1 
00  10  1=1. K 
TT(J)==T(N5) 

10  N5=N5*1 
L=<N*l>/2 
To=TT(D 

N6=N4 

N7=N4*1 

JU=1 

N2=N-1 

UN=!.C 

00  14  J=1.M2 
N5=N4-N1 

n?=n-j 

00  11  1=1. N3 
N8=N5*J 

TT(I!*{TT(I*1 j-TT(I) ) /(XT (N8)-XT (N5) ) 

11  NS=N5«1 

GOTO  (12.13).  JU 
It  UN=UN* fx-XT (N6) ) 

-•()=? 

N6=N6-1 
GOTO  14 

13  un=un*(x-xt:n7)) 

JU=1 
N7=N7  *  1 
L=L-1 
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TR=TR*UN»TT (L) 

fx=tp 

RETURN 

If (X-XT (1 J ) 16.16*17 
TR=FT  < 1 ) 

FX=TP 
RETURN 
TR=f T (NP) 

FX=TP 

return 

WRITE  £6.27)  NP ♦ NO 
GOTO  >5 

IF (X-2.*XT(1> ♦XT'?) )20. 15.15 

IF  (X-2.*XT  (NP)  *XT  (NP-1 )  US. 21.21 

2F(NP.LT«10)GOTO  23 

NS=NP~N 

N5=N5/2 

N!6=N^*N5 

If (XT(N6? ,GT.X)N4=N6 
IF (NS.GT.l)GOTO  22 
If  (X-XT(NA5)24»24.R 
If  (NA-N.'s!  25.9.25 
NA=NA»1 
GOTO  23 

WRITE  (6.28)  XT ( 1 ) 

GOTO  15 


FORMAT (21H  TABLE  TOO  SMALL  NP=  .IS.6H  ND=  . I5.2X .6HDVUINT) 
F0RMAT(23H  CONSTANT  TABLE  XT(1)=  .E14.7.2X.6HDVOINT) 

END 

SUBROUTINE  FHGKTiAX.AV.A7,VX,VY.VZ.XT,YT.ZT.OT) 

VXsVX.AX^DT 

VV=VY*AY«OT 

V7=VZ»AZ®0T 

XT=XT*YX*OT 

YT=YT*VY*OT 

ZT=Zr«VZ*F)T 

Rf TURN 

END 

SUBROUTINE  FlOfB.C.ZT.Z. VALUE. FIGS) 

DIMENSION  Z(6) .F1GM6.6) «C(6) 

8=ARS(B) 

1=0 

00  IS  N-1.5 

Ic  (B.C-F.C(N)  .AND.B.LT.C(N»i))  I=N 
CONTINUE 

IF(I.fO.0)GOTO  20 

GOTO  30 

CONTINUF 

GOTO  37 

J=0 

On  35  K=! ,5 

IF  (ZT  .GE.Z  (X!  .  AN0.2T  .LT.Z(N»1)5  J"*K 
CONTINUE 

IFU.EC.OIGOTO  20 

CALL  INTfP (FIGS.Z0.Z.8.C. J. I. VALUE) .RE TURNS (36 .37) 

Pf TURN 
VALUE=0. 

return 


SUBROUTINE  INtEP (FIG5«Z0tZ*3»C»J»I * VALUE?  .RETURNS (M»N) 

OIKENSION  FIG5(6»6).C(6J,Z(6) 

IF{FIG5<J*1»IO>.EQ.O.>GOTO  20 
IFJFIGSJ J*1 .1) .EO.O.JGOTO  20 
SH1*FIG5< J»I*1 ) 

SH2=FIGS<J,I) 

SH3*F!G5(J*1.!*1> 

SH4=FIG5(J*1 .1) 

Bi=C(I) 

S?=C(I*1> 

Z1=Z(J) 

Z2=Z  < J*1 ) 

3=ABS<£0 

SHft=SHl*{ <SH2-SM1)/(B2-01) >«<B-81) 

SHB=SH3« ( (SHA-GH3) / (82-R1 ) ) *  <8-81 ) 

VALUf=(ZO/(Z!*7?> )*fSHR-5HA>  *SH1 
VAt.UE-VALUE’1008. 

RETURN  N 
?0  CONTINUE 
VALUE=0. 

RETURN  m 
END 

SUBROUTINE  INTERP { XX ♦  I .  J.  A .N'.BFTA . EL » VALUE ,  YX i 
DIMENSION  A{5.17.6)«YX(17; ,EL(6) 

IF(A(N.J»1.I»} } .EO. 0. > GOTO  20 

AvF=A(N«J,I>»(JAfN,J*I»I)-AJN.J,I)  )  /  ( YX  i  J-M  )-YX  {  J)  )  )  «  t  XX-YX  (  J)  ) 
FYA=A<N,  J.l4l)*f(A{N,j4l,l*l)-A(N»J«I»l)  >/!YXiJ*l)~YX«J)  ))MXX-YX( 
1J)  5 

VALUE= ( (8ETA-EL  f I ) > / CFL  C I  * I ) -£L ( I ) > ) *  < A VF-FVA 5  *EVA 
RETURN 
20  CONTINUE 
VALUE-0. 

RETURN 

END 

SUBROUTINE  OROFR  (X> 

DI“ENSION  X (S) 

XK=1 
00  »C 
L=5-J 

DO  10  *=]  ,t 

!f{X(K)-*{K«miO.I0.6 
6  IFIXtKK) ,LE.X<KJ)KK=X 
10  CONTINUE 

DO  20  JJ=].S 
IFUJ.EO.KKJGOTO  20 
X«JJJ=X(JJ1®0.3 
?0  CONTINUE 
RETURN 
END 

SUBROUTINE  RELPOS (X»Y»Z*XT.YT*ZT » AX.BET A  »D»VX« VY# VZ) 

PJE*3, 1*1593 
XX=XT-X 
YY=YT-Y 
ZZ=ZT-Z 

D=SORT {XX*«*2*YY**2*ZZ**2) 

VV=SORT  (VX**2*VY**2»VZ««2) 

TT=SQRT  f  XT*«2->YT**2*ZT*»2) 

IF(VV. EO.O.. OR. TT, EO.O.JGOTO  30 
AX=ACD$< (VX»XT»VY«YT*VZ«ZT> /<VV*TT> ) 

SO  CONTINUE 

fi£TA=ASIN(ZZ/D! 

GOTO  15 
30  Ax=0. 
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,y#y-v*  *aWHU»  Win/* 


J" : 


^jsegrs? 


21  IF (D.NE.O.IGOTO  50 
3£TA=PIE/2. 

15  CONTINUE 
BFTA=BETA*57.?957 
AX=AX»57.2957 

IF(XT.rQ.0..AND.Yr.EO,O.)WPITF  (6.16) 

16  FORMAT (16H  TARGET  OVERHEAD) 

RETURN 

END 

SUBROUTINE  SHADOW { XX  « A. BETA .EL  .D*OELTAT .T.TO.EIGS.AO.PHI . W30.XT  .  YT 
1  *Z7*X0»Y0»Z0»YX«C) 

0 1 MENS ION  Z<6) .C(6> .EIG5 (6.6) .EL (6) »YX{1?) .4(5.17,6) 

DIMENSION  SIG (5* 1 0 »S) »  AuB (5) .A0{5) ,0ELAD0(5) .AH(5) .00(10) «A>.(5> 
COMHON/niMCOM/TABOIH  (120) 

SK= (T-70) /ALOO10 (ZT/ZO) 

8T=4.*SK6U./10o**4) 

»W=1 .5*W30*COS(PHI/57.?957) *(l./i0.**4) 

8=ST*BW 

CW=W30*COS(PWI/57.2957) 

IF(CW.GT.O.)OOTO  20 
IE<CW.£O.O..ANO.PT.LT.O.)GOTO  10 
IFfCW.LT.O.JGOTO  10 
GOTO  20 

NO  SHADOW  PEG ION 

10  CALL  ALPHO(XX.a.RETa.-L.O«OELTAT,T.TO.ZO.ZT.FIG5.AO.YX.Z.PHI.8.C.W 
630.xT*yt,xG.YO) 

RETURN 

?G  KCRW.EO.  (-BT)  )GOTO  21 
GOTO  22 

?1  xxx:=<-BT>/(l.S*W30*(l./lQ.**4S ) 

IF(XXX.GT.1.0.0».«XX.LT.-1.0)GOTO  10 

PHIC-ACOS(XXX) 

o.H  I  C=PH  I C  •  57  »  2957 

IF  (PHIC.LT. <50., AN0.PHIC.GT. 270.) GOTO  10 
IF ( (PHJ-PhIC) »LT.60.) GOTO  10 

22  CONTINUE 

IFCABS(R) .GT.CC6) .OR.ARStB) .LT  »C( i ) >GOTO  10 

CALL  AL (XX.A.SETA.FL.O,OELTAT,T.TO.ZO,ZT.FIG5-AO.YX,Z.PHI.8,C.W30, 
6XT.YT.X0.Y0) 

RFTU3N 

END 

FUNCTION  TA3 ( I . J.K.L  > 

C0MM0N/PIC0H/TA3DIW  (120) 

FUNCTIONAL  F0»  COMPUTING  MULTIPLY  OR  N  INDICES 
INOEX=J  *5-°  ( J-!  »4«  C<-1*3*  <L-I » )  ) 

tarstabdjmi index) 

I -INNER  MOST  fNDEX - 5  IS  ITS  NUMBER 

J-NEXT  INDEX - 4  IS  ITS  NUHREP 

K-NEXT  INDEX— 3  IS  ITS  NUMBER 
L-OUTERMOST  INDEX— IS  COMPUTED 
RETURN 

end 

SUBROUTINE  TABLE* (N.D.RETA.mBT.SBw) 

COMMON/D  I  MCOM/TAflD  I M  (120) 

IF1DETA.lt. 2. )J=0 

IF  (BETA. GE. 2.. AND. BETA. LT. 3.) J=1 

IF  (BETA.GE .3..SND.BETA.LT.5.) J=? 

IF (BETA. 6E. 5., AND. BETA. LT. 7. 5) J=2 
IF‘BETA.GE.T.S.ANO.BETA.LT.10.5) J=3 
IF <  J.FG.OGOTO  SO 
GOTO  55 
SO  CONTINUF 
8«T=1.0 


55 


10 


15 

20 

33 


RETURN 

1=0 

IF  (O.GE.l 000. .AND.D.LT .1500.5  1  =  1 
IF(O.GE.lS00..AND.D.t-T.2000»)I=2 
IF (D.GE.2000. . AND.D.LT .2500. > 1=2 
IF (O.GE. 2500.. AND.O.LT. 3000, >1=3 
IF (O.GE. 3000.. AND.D.LT. 3500.1 1=3 
IF 'D.GE. 3500.. AND. D.LE. 4000.)  1=4 
IF ( I .EQ.O) GOTO  50 
8BT=TAB (N» I  * J« 1 ) 

8BW=TA8 (N.I.J.2) 

I F (TAB CN. I «J»1> .EQ.1.0) GOTO  50 

RETURN 

END 

SUBROUTINE  VEG(8ETA,D.AT,N,AC,£«EE) 

DIMENSION  AT (5) • AC (5.4) ,E(4) ,EE(4) 

le (PETA.GT .£ (4) )GOTO  IS 

NW=0 

DO  10  1=1.3 

IF<8ETA.G£.£<I) .AND. BETA. LT.E ( 1*1 > )NN=T 
CONTINUE 

IF (NN.EO.O)GOTO  15 

C=  ( (BETA— £  (NN)  )  /  (E  (NN*  1 )  -E  (NN)  )  )«  (EE  (NN*I  )-EE  (NN) )  *EE  IW) 

GOTO  2C 

C=0. 

00  30  1=1.5 

AT ( I ) =AC ( I *N) »C* (D/1 000 . ) 

R-TURN 

END 

SUBROUTINE  WEEOIWANO, WIND. FEET. XT. YT.2T.X.Y.Z. PHI >W30) 
DATA  A/57. 2057/ 

W30=WIND/(3...?7*ALOG10(FEET/30.)) 

WX=SIN(WANO/AJ 

WY=CDS(WANO/A) 

Of>7=W/»  (X-XT)  ‘WY®  (Y-YT) 

B=SORT< (XT-X)o®2»(YT-Y>®®2« (2T-Z>®«2> 

PHI=AC0S(00T/8)*A 

RETURN 

END 
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